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We study the nonequilibrium energy transport across a topological insulator/superconductor junction, by 
deriving an interfacial heat current formula through scattering wave approach. Several anomalous thermal 
properties are uncovered, such as thermal energy's Klein tunneling, asymmetric Kapitza resistance and negative 
differential thermal resistance. We expect these findings could have potential applications for the energy control 
in various hybridized mesoscopic systems. 
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Topological insulators (TIs), characterized by a bulk gap 
and a gapless surface mode with a Dirac-Uke hnear dispersion, 
are presently one of the most interesting topics in condensed 
matter physics [1, 2]. Their conducting surface states in the 
insulating gap are topologically protected by time-reversal 
symmetry, hence robust to disorders and perturbations, po- 
tentially leading to various device applications. With the help 
of doping-induced superconductivity in TIs [3] or by deposit- 
ing superconducting materials on TIs due to the proximity ef- 
fect [4—6], the interplay between the superconducting ordering 
and the gapless chiral surface state has triggered much inter- 
est [2]. 

On the one hand, such topological insulator superconduc- 
tor (TI/S) junctions have been used to create chiral Majorana 
fermions for topological quantum computations and for the 
study of their impacts on electronic tunneling properties [7- 
12]. On the other hand, in view of the fact that understand- 
ing thermal properties of hybrid mesoscopic structures is cru- 
cial for future quantum/nano technology [13], hybrid topolog- 
ical superconductor junctions could be a natural candidate for 
managing heat dissipation/refrigeration and controlling en- 
ergy transport. It has the potential to open up a rich variety 
of thermal device concepts in superconducting-circuit-based 
systems, just like thermal diodes, thermal transistors, thermal 
logic gates and thermal memories in phononics [14]. 

In this Letter, we study the nonequilibrium energy transport 
across a TI/S junction interface and uncover its anomalous 
thermal properties, such as thermal energy's Klein tunneling, 
asymmetric Kapitza resistance and negative differential ther- 
mal conductance (NDTC). Among them, Kapitza resistance 
measures the interfacial thermal resistance when thermal en- 
ergy flows through the interface between two different mate- 
rials [15, 16]. Asymmetric Kapitza resistance is one unusual 
thermal property that the interface acts as a good thermal con- 
ductor if a positive thermal bias is applied, while with a nega- 
tive thermal bias it exhibits poor thermal conduction, thus ef- 
fectively acting as a thermal insulator. As such, it functions as 
a thermal rectifier or diode (for a review, see [14, 17]). While 
the NDTC, another unusual thermal transport phenomenon, 
where the heat current across a thermal conductor decreases 
when the temperature bias increases (for a review, see [14]), is 
an essential element for the construction of thermal transistors 
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FIG. 1. (color online). Schematic illustration of the reflections and 
transmissions at the TI/S interfaces. 



and thermal logic gates, and has been shown to exist in many 
nonlinear systems. We believe similar findings in supercon- 
ducting hybridized mesoscpic junctions could extend the con- 
ceptual thermal devices in Phononics into superconducting- 
circuit based systems and have great potential applications 
[18] in the near future. 

As depicted in Fig. 1, we consider a two dimensional (2D) 
TI/S junction attached to their respective thermal reservoirs at 
Tl and Tn, with an insulating barrier locating at x G [0,d]. 
The 2D TI junction part (x < 0) could be formed on the sur- 
face of 3D topological insulators, the S junction part {x > d) 
could be induced via the proximity effect, and the barrier (B) 
part could be constructed by gating voltage or doping. We as- 
sume the same Fermi level in both TI and S, and employ the 
Bogoliubov-de Gennes (BdG) equation = Etp to study 
the thermal transport properties. The Hamiltonian of the sur- 
face state on the topological superconductor is given by 



H = 



-A* 



A 



(1) 



acting on a Nambu basis {ipet, i>ei, il^ht, V'h), with -ffo(k) = 
hvp^k^dy — kyO^) + UQQ{x)Q{d — x) and A ^ 
iayA{9,T)Q{x — d). Here, vp is the Fermi velocity, (Jx{y) 
denote Pauli matrices, Uq is the barrier potential and A{6, T) 
depicts the order parameter with a given pairing symmetry and 
temperature dependencies. Throughout the work, we consider 
Ep ^ (|A|,iiJ) to satisfy the mean-field nature of BdG ap- 
proach, i.e., the superconducting coherence length hvp/\A\ 



2 



is much larger compared to the Fermi wavelength hvp/Ep- 
Adjusting Ep can be achieved via doping or gate voltage. 

Since in superconductor junctions, the quasiparticles are 
carriers of thermal energy, we need to obtain the transmis- 
sion (equivalently the reflection) coefficients of quasiparti- 
cles in order to investigate the thermal transport properties. 
Considering the conservation law for particle current, we can 
simplify the problem by merely considering the particle (in- 
cluding both electron and hole) current in the side of Tl. 
Defining / — (^"p and its hole counterpart g — (^J^p, 
we have the probability density for finding either an elec- 
tron or a hole, P — |/p + jgp. By using the BdG equation 
ihdfjjg) — -ff (g) with Eq. (1), and considering the continuity 
equation dtP + W ■ Jp = 0, we obtain the x-component of 
the particle current: 

Jp = vpiP&yf - g^'&yg) . (2) 

Note the hole current is a time-reversed counterpart of the 
electron contribution so that it naturally obtains an oppo- 
site sign compared with the electron current. If we express 
the whole wave function at the Tl side in terms of the nor- 
mal reflection amplitude, b, and Andreev reflection [19] (for 
electron-hole conversion) amplitude a, which will be defined 
explicitly below, the particle current can be written as 

Jf,{E,9) = 2vfcos9{1 - \a{E,e)\^ - \biE,0)\^) ■ (3) 

This expression of the particle current has a clear physics pic- 
ture that fi{E,0) := 1~ \a{E,9)\'^ - HE,0)\'^ denotes the 
transmission of energy carriers with energy E and incident 
angle 9, vp cos 9 the effective velocity in the x direction and 
2 the spin degeneracy. Considering the carrier's energy E, the 
Fermi occupation difference between two sides of the inter- 
face fL-fR = ^eiuIt^^^ - ^Eik^Tj^^^ and summation over 
all possible incidence angles and momenta /"^^j '^^ ^ 

\dEfJ^\d9{dEldK)-^ = JdEj:i%d9j-^, we 
arrive at the energy current expression: 

2 f°° 

Jq = IJ dEEK{E)[fL ~ fp] , (4) 

with k(£;) = J^l^^^d9 COS 9k{E, 9) = J^(^^^^d9 cos 9(1 - 
I a p — 1 6 P ) . Note that this expression is general and is obtained 
before solving the reflection coefficients a,b. In fact, similar 
expressions of the energy current in topological-trivial metal 
superconductor junctions have been obtained by the rigorous 
derivation through linear response of entropy production [20] 
or by heuristic argument [21]. It is worth emphasizing that 
all the transports considered in this work are charge neutral, 
i.e., the carriers transport only thermal energy without charge 
current. In fact, if we follow the same procedure for charge 
transport, we will arrive at a similar expression for electric 
current: J, = ^ J dE J d9 cos 9a{E,0)[fL - fp], with = 
1 + lap - similar to the BTK formula [19, 22]. One 
can then get that since the even symmetry of a {a{E,9) = 



a{—E, 9)) and the odd symmetry of fp — fp with respect to 
E, Je vanishes as zero. 

We now proceed to determine the scattering coefficients 
of Andreev reflection amplitude a and normal reflection am- 
plitude b by imposing the boundary conditions on the wave 
functions at the interfaces of barrier. Diagonalizing Eq. (1) 
straightforwardly yields the wave functions in the Tl, barrier, 
S regions. In the Tl region {x < 0), for electrons and holes 
traveling the zLx direction with a conserved transverse mo- 
mentum ky and energy E measured from Ep, the wave func- 
tions are given as 

where k^ — Epcos9/{hvF) and 9 = aicsm{hvpky/Ep) 
is the angle of incidence. Note we have used the condition 
Ep ^ E. In the barrier region (0 < x < d), we employ the 
thin barrier limit [23]: d 0, Uq —> oo but with a constant 
product Z = Uod/{hvp), characterizing the strength of the 
insulating barrier, we then obtain 

V;^'± = (0,0,i,±l)e*(=F^^/'i+fe«w). 

In the S region {x > d), the electron-(hole-)like quasiparti- 
cles are mixtures of electrons and holes. Thus, the transmitted 
wave functions have the forms: 

where e"^± = A{9±,T)/\A{9±,T)\ with 9+ =9,9^ = 
TT — 9 and r± = v±/u±, with u^. — ^(1 + 
^E^ - |A(6'±,T)|V|£;|) = 1 - w|. For a d-wave pair- 
ing symmetry, A{9,T) = A(r)cos(26l - 2a) with A(T) 
Ao tanh[(7rfcBTc/Ao) V0-953(7;/T - 1)] [24]. Here is 
the critical temperature, Ao denotes the superconducting gap 
at zero temperature and a is the angle between the normal 
direction of the barrier interface and the x axis of the d^2_y2- 
wave superconductor. By taking into account the boundary 
conditions: 

^'tiU=0 = *bU=0, *BU=d = *sU=d, 

with *Ti = i^jt + bi^ji + aV^^r, *s = ieV's^ + thi^s'^, 
vI^B = ri7/iB+ + r2i/Js^ + J'sV'b^ + f47/'B , *e Andreev and 
normal reflection coefficients are found to be 

" ~ cos2 9 + sin^ Z sin^ 9{1 - r+r_e»('*--'^+)) ' 

^ _ sinZsin6'(cosZcos6' - isinZ)(l - r+r_e*('^--'^+)) 

~ -e-*«[cos2 9 + sin^ Zsin^ 61(1 - r+r_e»('*--'^+))] 
Finally, using the obtained coefficients a and b, we get 
ii{E,9) : = 1 - \a{E,9)\^ - \b{E,9)\^ 

i-«Ti|r+p + (/^Ti-i)|r+r_p 
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(b) c(/(it/4) 
FIG. 2. (color online). Thermal energy's Klein tunneling, (a) Ther- 
mal conductance as a function of the barrier strength Z and the rota- 
tion angle a of the d-wave order parameter, (b) Intersections of (a) 
for Z — 0, and 7r/2. (c) Angle-resolved thermal conductance Gg for 
different a and Z. Solid lines are for Z — while dashed ones for 
Z = 7r/2. Parameters are T = 70K, Tc = 116K, Ao = 30meV. 
Thermal conductance G is in the unit of Aq/cb /h. 



where 



KTi:-^(|-B|»Ao,0) = 



cos^e 



sin Z sin" 9 



(6) 



is reminiscent of the relativistic Klein tunnehng [25], as a con- 
sequence of the spin-orbit coupHng in TIs. The barrier be- 
comes transparent for the thermal energy transport at the res- 
onance condition, Z :— — m:,n — 0,±1, ■■ ■ (such that 
sin Z — 0), or at the normal incidence (9 — 0). 

In the linear response regime ^ T + ST/2,Tjj^ = T - 
ST/2, 5T 0, we have the thermal conductance: 



G 



Jq 

ST 



dE 



-n/2 



cos 0k{E,0) 

V2 ' 4fcBT2 C0sh2(2j|^ 



d9- 



(7) 

As illustrated in Fig. 2(a), the oscillatory dependence of G on 
the barrier strength Z shows the Klein tunneling behavior of 
energy: thermal conductance anomalously increases when the 
barrier increases within [(n+l/2)7r, (n+l)7r], consistent with 
the behavior of sin^ Z in Eq. (6). G also has an oscillatory 
dependence on a, the rotation angle of the superconducting 
order parameter. As detailed in Fig. 2(b), when Z = 0, G 
monotonically decreases as a rotates from to 7r/4, but when 
Z increases to 7r/2, G first increases and then decreases. This 
nonmonotonic behavior is due to the competition of Klein tun- 
neling and superconducting order parameter (gap), which we 
explain below. 

In the S pait, thermal energy is carried by the quasipar- 
ticles, which only transport beyond the gap A, so that the 
smaller A, the larger thermal conductance. In other words, 
the rotation angle a adjusts G through adjusting the angle de- 
pendence of the superconducting gap. To further understand 
the competition between Klein tunneling and order parameter 
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FIG. 3. (color online). Rectification ratios as a function of tempera- 
tures. Values apart from 1 indicate asymmetric Kapitza resistances. 
■ (T = 35K), • (T = 70K), A (T = 104K), T (T = 140K) denotes 
for Z = 0. Their hollow counterparts are for cases of Z = 7r/2. In- 
set shows one example of the Jq profiles depending on temperature 
bias for Z = (soUd) and 7r/2 (dashed), at T = 70K. Other param- 
eters are the same as in Fig. 2. The energy current Jq is in the unit 
ofAl/h. 



A, we define the angle-resolved thermal conductance Ge as 



G 



— /^^/2 d9Gg and plot it as a function of the incidence 
angle 9 in Figs. 2(cl-c3). For the case of Z = 0, when 
a = 0, Ge has two peaks around 6 = ±7r/4 where A is 
gapless; while, when a rotates to 7r/4, two peaks becomes a 
single peak around 6 = 0. As a result, G, the angle inte- 
gration of Gg, decreases as a increases from to 7r/4. For 
a changing from tt/A to 7r/2, the behavior is symmetrically 
reversed. When the barrier Z increases, the Klein tunneling 
comes into play. For the case of Z = tt/2, when a = 0, 
although there are still two peaks for thermal conductance, 
their intensity is suppressed dramatically. When a = tt/A, the 
gapless angle coincides with the normal incidence angle, the 
barrier becomes transparent. Therefore, the Klein tunneling 
helps to keep the single conductance peak (at 9 = 0) intensity 
unchanged. At the intermediate regime a = vr/S, one peak 
near 6 = preserves while the other peak far from 9 = 
is repressed. As a consequence, the thermal conductance G 
increases first and then decreases within a £ [0, 7r/4]. 

The asymmetric Kapitza resistance is essentially a nonlin- 
ear response behavior. It is a consequence of different tem- 
perature responses of different materials at two sides of the 
interface [14]. The inset of Fig. 3 shows typical Jq behaviors 
at T = 70K for Z = (soHd line) and 7r/2 (dashed line), 
via varying temperature bias. As a measure of the asymme- 
try, we define the rectification ratio R = | J_/J+|, where J^ 
refers to thermal current when ST = Tl — Tu > while 
J_ refers to thermal cuiTent after switching the temperature 
bias Tl ^ Tff . As shown in Fig. 3, when temperature bias 
is apart from 0, R becomes more deviating from 1. Except 
for the low temperature case, (e.g. T = 35K), other three 
examples show that increasing bias does not always increas- 
ing the rectification ratio at the large bias regime. The results 
also indicate that although the insulating barrier changes the 
Jq profiles quite noticeably, it does not change the rectifica- 
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FIG. 4. (color online). NDTC for various temperatures at Z — 
(filled symbols) and Z — ■k/2 (hollow symbols). Other parameters 
are the same as in Fig. 2. Results imply that the barrier can enhance 
the NDTC. 



tion ratio significantly, which is even slightly enhanced by the 
barrier, hi addition, R has a nonmonotonic temperature de- 
pendence that increasing T first increases and then decreases 
R, as exemplified by the highest curves for T ^ 70K in Fig. 
3. This is reasonable that the asymmetric Kapitza resistance 
in our system results from the different temperature responses 
of topological insulator and superconductor at the sides of the 
interface: At higher temperature, the superconducting gap di- 
minishes so that the superconductor tends to topological in- 
sulator and both sides of the interface tend to have the same 
temperature response, which explains the reduction in R. At 
lower temperatures, although two sides of the interface have 
distinct temperature responses, the bias ST cannot be larger 
than T such that small bias reduces R. Therefore, the optimal 
R appears at intermediate temperature. 

The superconductor is usually a bad thermal conductor 
since its gap A forbids the existence of quasiparticles that is 
responsible for the energy transport. Increasing the tempera- 
ture Tji at the S part could increase the energy transmission 
k{E) (in Eq. (4)) by diminishing the superconducting gap to 
allow more thermal energy carriers, while simultaneously de- 
crease the temperature bias ST = T^—Tr as well as the occu- 
pation difference Jl- fR= ^e/u^t^^^ - -^WTkh^v ^"'^^ 
the increased k(-E) is able to compensate the loss in — //f 
due to the decreased ST, we expect to observe NDTC, /.e., the 
energy current increases as the bias ST decreases. To char- 
acterize this anomalous behavior, we define the differential 
thermal conductance (DTC): dJqjdST, for given finite ST. 
Note that this is a nonlinear quantity similar to the nonlin- 
ear differential electrical conductance dl /dV, and is different 
from the previous defined thermal conductance G in Eq. (7), 
which is a linear quantity at ST — >■ 0. As a showcase, we 
calculate DTC in Fig. 4, with varying Tr in the S region but 
fixing Tl as a reference temperature. Indeed, NDTC appears 
as we expect, although it is absent for the low Tl case {e.g. 
Tl = 70K > Tr). Increasing enhances the regime of 
NDTC. However, when Tr increases across a threshold and 
approaches to Tl, NDTC disappears. In addition, the insu- 



lating barrier can also slightly enhance NDTC, as shown in 
Fig. 4. When the superconducting part is replaced with a topo- 
logical insulator, there exists no NDTC. 

In summary, using the scattering wave approach, we have 
deiived an interfacial heat cuiTent formula in a TI/S junction. 
With the help of this formula, we have studied the nonequilib- 
rium energy transport across this interfacial system and have 
uncovered several anomalous thermal properties for the TI/S 
interface, such as thermal energy's Klein tunneling, asymmet- 
ric Kapitza resistance and negative differential thermal resis- 
tance. We believe our results have revealed the potential ap- 
plications of hybridized superconductor/topological insulator 
systems for the smart energy control at mesoscopic scales. 
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